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Introduction: Clay minerals have been identified 
in the central peaks and ejecta blankets of impact cra- 
ters on Mars [e.g. 1,2]. Several studies have suggested 
these clay minerals formed as a result of impact- 
induced hydrothermalism either during Mars’ Noachi- 
an era or more recently by the melting of subsurface 
ice [3-6]. Examples of post- impact clay formation is 
found in several locations on Earth such as the Mjolnir 
[7] and Woodleigh Impact Structures [8]. Additionally, 
a recent study has suggested the clay minerals ob- 
served on Ceres are the result of impact-induced hy- 
drothermal processes [9]. Such processes may have 
occurred on Mars, possibly during the Noachian. 

Distinguishing between clay minerals formed pre- 
or post-impact can be accomplished by studying their 
IR spectra [10,11]. In fact, [10] showed that the IR 
spectra of clay minerals is greatly affected at longer 
wavelengths (i.e. mid-IR, 5-25 pm) by impact-induced 
shock deformation while the near-IR spectra (1.0-2. 5 
pm) remains relatively unchanged. This explains the 
discrepancy between NIR and MIR observations of 
clay minerals in martian impact craters noted in [12]. 
Thus, it allows us to determine whether a clay mineral 
formed from impact-induced hydrothermalism or were 
pre-existing and were altered by the impact [10]. Here 
we study the role of impacts on the formation and dis- 
tribution of clay minerals on Mars via a fully 3-D 
Monte Carlo cratering model [13,14], including im- 
pact-melt production using results from modern hy- 
drocode simulations [15]. We identify regions that are 
conducive to clay formation and the location of clay 
minerals post-bombardment. 

Impact Model: Mars is modeled as a Cartesian 
sphere of radius 3390 km discretized into cubic vol- 
ume elements 10 km on a side. We use Monte Carlo 
methods to select a population of impactors with a 
size-frequency distribution similar to the present-day 
asteroid belt and model the bombardment following 
the E-belt hypothesis [16,17]. The total bombardment 
mass incident on Mars is constrained using the esti- 
mated lunar bombardment mass, -3.5 xlO 19 kg [17], 
assuming Mars receives -2.76x more material than the 
Moon [18]. Though it receives more impactors due to 
its proximity to the asteroid belt, the mean impact ve- 
locities on Mars is -0.54* that of the Moon [18]. Im- 
pact velocity is thus modeled following a Rayleigh 
distribution about a mean value of 7 km/s and 1 1 km/s 
for the pre- and post-LHB eras, respectively. 

Once chosen per run, the synthetic impactor popu- 
lations are converted to crater populations following a 


Pi scaling law [19], where a projectile of diameter D 
and density pi impacting at velocity v/ produces a 
crater of transient diameter 
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where (pi/p m ) - 1, g = 3.7 m s" 2 is Mars’ gravity, and Q 
is impact angle, which is randomly selected following 
a distribution of dQ=sin2Q. Impact location for every 
projectile onto the simulated Mars is chosen randomly 
for longitude and for latitude follows d4>=sin24>. An 
impactor excavates and ejects material within a vol- 
ume approximated as an oblate spheroid of depth 
D t J 8, which is emplaced within an ejecta blanket 
modeled as an annulus surrounding the transient crater 
extending D tc from the impact point. Impact melt is 
modeled following the scaling relationships from [15] 
such that an impactor of radius r p produces a volume 
of melt at a depth of 
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where a£ = 2.52 and b^ = 0.651, of radius 
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where a^ = 3.00 and b^ = 0.674. Within the impact 
melt, clay mineral formation can occur and over time 
be redistributed across the surface. 

Preliminary Model Results: For impactor popu- 
lations with SFD characteristics like the asteroids belt, 
the bulk of the impacting mass is characterized by 
D-100 km, of which there are -70 such impacts on 
Mars. Following Equations 2-3, the bulk of the impact 
melt volume would then be characterized by z-50 km 
and r-55 km such that a total of -5><10 7 km 3 impact 
melt volume is produced during Mars’ bombardment 
history. This is equivalent to -0.2 km think global lay- 
er of melt. In Figure 1, we show an example simula- 
tion demonstrating the post-bombardment distribution 
of impact-melt on Mars, where P me it is the volumetric 
probability of finding impact-melt. By our assump- 
tions, this is the most likely area to find impact-formed 
clays. Though large impacts produce the most melt, 
their large excavation volume implies efficient mixing 
with undisturbed material. Smaller impacts, though, 
can more easily eject these melt volumes [20]. 
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Figure 1: Post-bombardment distribution of impact-melt on Mars from a single impact history. The colors represent Pmelt. 
Crater floors, where mass-wasting, which is not included in the model, likely modifies the composition, are denoted by black. 


For this simulation, on average P me it = 0.03. Thus, 
though Mars receives a lot more impacting mass than 
the Moon, the small impact velocity combined with the 
large surface area result in a small probability of find- 
ing impact melt relative to the Moon. 

Implications for Mars: Clay minerals have been 
identified by their NIR spectral signature in the central 
peaks and rims of impact craters on Mars [e.g. 1,2,21]; 
however, this spectral signature results from only the 
first few pm depth of the surface and the absolute vol- 
ume of the clay mineral is not certain. In our model, 
where there is impact melt, there is a possibility of 
forming clay minerals, but not all of that model vol- 
ume will result in clay mineral formation. Our mod- 
eled P me it is thus an upper bound probability of finding 
clay minerals. 

Though most melt is made by large impacts, a large 
amount of undisturbed material is also excavated, so 
the ejecta is a mixture of both. Smaller, later impacts 
near these large basins excavate the basin melt and 
because their excavation volume is small is dominated 
by it resulting in high P me it (red areas in Fig. 1). The 
resulting higher P me i t around smaller, younger impacts 
suggests areas of higher probability of clay mineral 
formation. In fact, a comprehensive study by [22] 
showed that given its slightly alkaline conditions (pH 
6-8), hydrothermal systems resulting from terrestrial 
impacts are more likely to form clay minerals in small- 
er craters (D < ~65 km). 

Future Work: The model result presented here are 
only one possible scenario in which clay minerals can 


result from post- impact hydrothermal activity. We 
expect that different starting conditions (e.g. target 
rock composition, water content) in the model will 
yield different results. From this further work, we ex- 
pect to be able to constrain the conditions under which 
clay minerals can form from impact-induced hydro- 
thermal systems and the (re-)distribution of such clay 
minerals in the impact ejecta. 
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